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Context, The spectral energy distribution (SED) in chemically peculiar stars may be significantly affected by their abundance anoma- 
lies. The observed SED variations are usually assumed to be a result of inhomogeneous surface distribution of chemical elements, 
flux redistribution and stellar rotation. However, the direct evidence for this is still only scarce. 

Aims. We aim to identify the processes that determine the SED and its variability in the UV and visual spectral domains of the helium- 
weak star CU Vir. 

Methods. We used the TLUSTY model atmospheres calculated for the appropriate surface chemical composition to obtain the emer- 
gent flux and predict the rotationally modulated flux variability of the star. 

Results. We show that most of the light variations in the vby filters of the Stromgren photometric system are a result of the uneven 
surface distribution of silicon, chromium, and iron. Our models are only able to explain a part of the variability in the u filter, however. 
The observed UV flux distribution is very well reproduced, and the models are able to explain most of the observed features in the 
UV light curve, except for the region 2000 - 2500 A, where the amplitude of the observed light variations is higher than predicted. 
The variability observed in the visible is merely a faint gleam of that in the UV. While the amplitude of the light curves reaches only 
several hundredths of magnitude in the visual domain, it reaches about 1 mag in the UV. 

Conclusions. The visual and UV light variability of CU Vir is caused by the flux redistribution from the far UV to near UV and 
visible regions, inhomogeneous distribution of the elements and stellar rotation. Bound-free transitions of silicon and bound-bound 
transitions of iron and chromium contribute the most to the flux redistribution. This mechanism can explain most of the rotationally 
modulated light variations in the filters centred on the Paschen continuum and on the UV continuum of the star CU Vir. However, 
\q ' another mechanism(s) has to be invoked to fully explain the observed light variations in the u filter and in the region 2000 - 2500 A. 
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,— ' \ 1. Introduction tribution induced by these lines has been suspected to signifi- 

cantly affect the spectral energy distribution (SE D) and to be the 
Chemically peculiar (CP) stars are among the most enigmatic prime SQUrce of ±e Hght variability (e ^ Molnar 1973; Kodaira 

7^ . obj ects of the upper pail of the main-sequence. The processes of m) Bound . free m *ions (Peterson 1970 HLanz etal.l ll996) 

radiative diffusion and gravitational settling in the atmospheres were alsQ exp ected to play some role. Surface temperature dif- 

of these stars cause pronounced deviations from the solar value fences o r variable temperature gradients (Weiss et al. 1976; 

in the chemical composition dVauclaiij | 2003i | Michaud| | 2004|) . ^^T 9 78) were on the list of possible causes of the SED vari- 

g , Apart from the chemical peculiarity, many of CP stars show van- abm as wdl The preS e nC e f a magnetic field may affect the 

ations in the magnetic field as well as in their spectra and light. S E D and its v ariability, provided the field is sufficiently strong 

These variations are usually strictly periodical and modulated (Kochukhov et aL 2005). Finally, hot CP stars may have circum- 

by the rotation of the star. The uneven distribution of the sur- steUar sheUs fed by t he stellar wind, cau sing variability through 

face magnetic field (with dipole component dominating in most the Hght absorption (Landstree t & Borra 1978; Nakaiima 1985; 

cases) is one of the factors that cause the uneven distribution of |^ ith & Groot£ 2m . Xownsend et al. 20jlf 
chemical elements and, consequently, also the periodic spectrum 

variability. The uneven distribution of chemical elements is sus- Realistic SED simulations of CP stars were not possible un- 

pected to be the origin of the light variability, but this connection til the techniques of Doppler mapping and model atmosphere 

is still not very well understood. calculations were considerably developed. The Doppler map- 

The line blanketing caused by numerous lines of overabun- ping enables one to precisely map the distribution of individ- 

dant elements (mainly the iron-peak ones) and the flux redis- ual e le ments on the su r face o f rotating stars (e .g.. iRice et all 

1 19891; iKh okhlova et aLlBOOOt lPiskunov & Kochukhovl 120021: 

Send offprint requests to: J. Krticka, Liiftinger et al.l l2010t iBohlender et al.1 |2Q1Q|) . Detailed model 
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e-mail: krticka@physics.muni.cz atmospheres (e.g.. lLanz & Hubenvl 12007) enable one to pre 
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Table 1. CU Vir parameters from spectroscopy (Kuschnig et al. 
1999). 



Effective temperature 13 000 K 

Surface gravity log g (cgs) 4.0 
Inclination i 30° 
Rotational velocity projection v rot sin i 160 km s -1 

Helium abundance -3.2 < £He < -1 

Silicon abundance -4.6 < £si < ~2.3 

Chromium abundance -6.6 < £q < -4.4 

Iron abundance -5.5<SFe<-3.5 



cisely predict the fluxes from surface elements with a peculiar 
chemical composit ion, taking into ac count realistic bound-free 
dSeaton et al.ll9 92) and bound-bound ( Kuru czl 1 9941) transitions. 

Based on precise model atmospheres and detailed sur- 
face maps it was possible to follow earl i er attempts of light 
curve modelling dKrivosheina et al.l Il980b [Rvabchikova 1990) 
and t o successfully simu late the light curves of several CP 
stars. IKrticka et all (120071) showed that the light variations of 
HD 37776 (7^ = 22 000 K) can be explained to be a result 



of inhomogeneous surface distribution of helium and silicon. 
IKrticka etalJ d2009) showed that most of the observed light vari- 
ations in HR 7224 {T e g = 14 500 K) are caused by inhomo- 
geneous surface distribution of silicon and iron. F or the cooler 
star e UMa (r eff = 9 000 K) IShulvak et all (1201 Oal) showed that 
chromium can also contribute significantly to the light variabil- 
ity. Moreover, modern atmosp here models are able t o explain the 
observed SED in detail (e.g.. lShulvak et ai]l2010bl) . pointing to 
the importance of rare-earth elements in cooler CP stars. 

For our present study we selected one of the most enig- 
matic CP stars, CU Vir (HR 5313, HD 124224). The light vari- 
abilit y of CU Vir has been known for more than half a cen- 
tury (lHardidll956l) . CU Vi r belongs to a rare group of CP stars 
that show period changes dPvper et alJll998UPvper & Ade lman 
12004 iTrigilio et alJl2008ll201ll:lMikulasek et al.rt201lh . theoret- 
ically studied by Stepien ( 1998). Moreover, CU Vir is a source of 
variable radio emission, resembling a ra dio lighthouse of pulsars 
dTrigilio et alfcOOOLlKellett et all2007l) . 

All these observations make CU Vir one of the most appeal- 
ing targets for theoretical studies. We studied the nature of the 
SEP variations of thi s star using Doppler abundance maps of 
Kuschnig etall dl999l) . 



2. Simulation of the SED variability 

2.1. Stellar parameters 

The stellar parameters of CU Vir and abundan ce maps of he- 
lium, s ilicon, chromium, and iron adopted from Kuschnig et al.l 
(1999) are given in Table [TJ Note that there i s also a magnesium 
abundance map available in lKuschnig et al.l (1 19991) . but because 
of the low maximum magnesium abundance derived we did not 
include its inhomogeneous surface distribution. The abundances 
in the maps are expressed as log(N e \/N tot ), but we used abun- 
dances relative to hydrogen, i.e., e e i = log (N e i/Nn)- 

The calculation of the rotational phases for individual obser- 
vations is not a straightforward task, because the instant rota- 
tional period P(t) of the surface layers of CU Vir is changing. 
Accor dingly, we applied the new ephemeris of Mikulasek et al.l 
d201 ll) . where the phase function i?(f) is approximated by the 



fourth-order polynomial of time: 

A /, . .x t - Mq t -To 

m = #<>- — (\e 2 - & 4 ) ■ § Q = —A e = — °, (i) 

P(t)= l/&± P o [l+A/n(3&-40 3 )], (2) 



where P(t) is the instant period at the time f, i?o is the phase 
function for a linear ephemeris with the origin at Mp = 
2 446 730.4447 and the basic period P n. lMikulasek et al.l (1201 ll) 
found that P = 0^52069415(8), A = $5643(29), 77 = 
1 3260(70) d, and T = 2446 636(24). The formula accounts for 
the period variability observed in CU Vir and enables us to de- 
termine the rotational phase with an accuracy better than 0.002 
P. 

The phase shift between this ep hemeris </> = frac(ff(Q) (frac- 
tional part of #(f)) and that used by Kusch nig et al.l £1999) 0kus 
in the time of their spectral observations was A(p = <p - 0Kus = 
0.52045, and the origin of the phase function was determined to 
be at HJD 2446730.4447. 



2.2. Model atmospheres and synthetic spectra 

We used th e code TLUSTY for the mod el atm o sphere 
calculations dHubenvl [l988l iHubenv & Lanzl fl99l fl995t 
lLanz & Hubenv 2003 ). Although the code enables us to cal- 
culate NLTE models, we confined ourselves to the LTE plane- 
parallel models, because we expected the NLTE effects to be 
marginal for the light variability. The atomic data (taken from 
lLanz & Hubenvl 120071) were selected to be appropriate for B 
type stars, the atomic d ata for silicon in pa rticular are based on 
iMendoza et all ( [1991 . iButler et al.l (fl993b and Taylor ( 2011) 
in preparatiorfl for irononlKuruczldl 9941). iNaharldl 9961) . lNaha3 
dl997l) . lBautista & Pradhanldl997h , andlBautistal(ll996l) and fo r 
other elements onlLuo & Pradhanl d 19891). iFernlev et aL ( 1999 1 ) . 
iTullv et al.1 (fl990h. iPeach et al.1 (1 19881) . Iffibbert & Scottl i 19941) . 
and lNahar & Pr adhan (119931) . We prepared our own ionic mod- 
els for chromium (Cru-Crvi) using data taken from Kurucz 
(2009fl 

We assumed fixed values of the effective temperature and 
surface gravity (according to Table [TJ and adopted a generic 
value of the microturbulent velocity v tur b = 2kms~ 1 . The abun- 
dance of helium, silicon, chromium, and iron differed in individ- 
ual models as e xplained below. We used the solar abundance of 
other elements ( As plund et al.ll2005h . 

For the calculation of synthetic spectra we used the 
SYNSPEC code. The synthetic spectra were calculated for the 
same parameters (effective temperature, surface gravity, and 
chemical composition) as the model atmospheres. We also took 
into account the same transitions as for the model atmosphere 
calculations. To this end we included the same chromium and 
iron lines as we used for the model atmosphere calculation in 
our SYNSPEC line list. This is not particularly important in the 
visible, but in the ultraviolet (UV) numerous lines for which 
only the theoretical data are available significantly influence the 
spectral energy distribution. Additionally, we included the lines 
of all elements with the atomic number Z < 30, that were not 
accounted fot the model atmosphere calculation. We computed 
angle-dependent intensities for 20 equidistantly spaced values of 
fi = cos 8, where 8 is the angle between the normal to the surface 
and the line of sight. 



Note that the opacity caused by autoionisation is included via 
bound-free cross section (as default in TLUSTY). 
2 http://kurucz.harvard.edu 
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Table 2. Individual abundances Eu e , ss\, Sa, and epe of the 

ne, Lr ' re 20000 

model grid 



He 


-2.0 


-1.0 










Si 


-4.75 


-4.25 


-3.75 


-3.25 


-2.75 -2.25 


16000 


Cr 


-6.4 


-5.9 


-5.4 


-4.9 


-4.4 




Fe 


-5.4 


-4.9 


-4.4 


-3.9 


-3.4 





The model atmospheres and the angle-dependent intensi- 
ties I(A, 9, sue, ssi> sqi, £Fe) mentioned above were calculated for 
a four-parametric grid of helium, silicon, chromium, and iron 
abundances (see Table |2j. For silicon this gri d fully covers 
the ra nge of silicon abundances in the map of iKuschnig et alj 
(Il999h . but for helium, chromium, and iron the lowest abun- 
dances detected by Kuschnig et al.l ( 1 1999b are omitted from the 
grid. Our test showed that this restriction of the grid does not 
influence the predicted light curves. The generation of the com- 
plete grid would require the calculation of 300 model atmo- 
spheres and synthetic spectra. Because not all a bundance com- 
binatio ns are required for the interpolation of the Kuschnig et al.1 
( 1999) maps, we calculated only those models that were neces- 
sary. This helped us to reduce the number of calculated models 
by half. 



2.3. Phase-dependent flux distribution 

The radiative flu x in a colour c at the distance D from the star 
with radius R t is ( Mihalal [l978l) 



fc 



(1)7 



7 c (6»,Q)cos6>dQ, 



(3) 



visible 
surface 



where the intensity I c (6, Q) at each surface point with spherical 
coordinates Q is obtained by means of interpolation between the 
intensities I c (0, sue, £si> £&> epe ) calculated from the grid of syn- 
thetic spectra (see Table |2]l as 



I C {0, SHe, SSi, £Cr, £Fe) 



Jo 



(A) I(A, 6, £ H e, £Si, £Cr, S Fe ) &A. 



(4) 

The transmissivity function <S> C (A) of a given filter c of the 
Stromgren photometr ic system is approx imated for simplicity by 
a Gauss function (see lKrticka et al.|[2009l for details). 
The magnitude difference is defined as 



Am c = -2.5 1og(^), 



k 

IT 



(5) 



where f c is calculated from Eq. [3] and /J ef is the reference flux 
obtained under the condition that the mean magnitude difference 
over the rotational period is zero. 

3. Influence of the abundance on the emergent flux 

Individual elements modify the temperature distribution of 
model atmospheres by their bound-free and bound-bound transi- 
tions. This can be seen in Fig.[TJ where we compare the temper- 
ature distribution of model atmospheres for typical abundances 
found on the surface of CU Vir. The bound-free (caused by ion- 
isation of helium and silicon) and bound-bound transitions (line 



12000 ■ 



< 




Fig. 1. Upper plot: The dependence of temperature on the 
Rosseland optical depth r ross in the reference model atmosphere 
with e He = -1-0, esi = -3.75, sq = -5.9, and e Fe = -4.4. 
Lower plot: The temperature in the model atmospheres with 
modified abundance of individual elements minus the temper- 
ature in the reference model atmosphere. 



transition of chromium and iron) absorb the stellar radiation, 
consequently the temperature in the continuum-forming region 
(i"ross ~ 0.1 - 1) increases with increasing abundance of these 
elements. For silicon and iron the influence of abundance on the 
temperature is stronger, for chromium the influence is weaker, 
while for typical helium abundances found on the surface of 
CU Vir the changes of temperature are only marginal. 

In atmospheres with overabundant helium, silicon, 
chromium or iron the enhanced opacity leads to the redis- 
tribution of the flux from the short-wavelength part of the 
spectrum to the longer wavelengths of the UV spectrum, and 
also to the visible spectral regions (see Fig . [2j» . Consequently, the 
overabundant spots are bright in the uvby c olours, and are dark i n 
far-ultraviolet bands. As already found by Krtic ka et al.l d2007l) . 
helium can affect the flux distribution only if it significantly 
dominates over hydrogen, i.e. for Eu e > 0.5. Consequently, 
for model atmospheres with underabundant helium the flux 
variations are only marginal. Note also that the flux variations 
caused by silicon are most pronounced in the far-UV region 
with A < 1600 A. 

These flux changes can be detected as a change in the appar- 
ent magnitude. To demonstrate this, we plot (Fig. [3) the relative 
magnitude difference defined as 



Am A 



-2.5 log 



Hf 



(6) 



against wavelength. Here Hf { is the reference flux calculated for 
slightly overabundant chemical composition (with eee = -1-0, 
esi = -3.75, ecr = -5.9, and epe = -4.4). As can be seen 
in Fig. [3] the absolute value of the relative magnitude differ- 



3 



J. Krticka et al.: Modelling of the ultraviolet and visual SED variability in hot magnetic Ap star CU Vir 



7x10 
„ 6x10 
■< 5x10 7 



cj 4x10 
o 3x10 
& 2x10 7 

1 x -\0' 



< 




1000 



2000 



3000 



4000 



5000 



MA] 



6000 



Fig. 2. Upper plot: The emergent flux from a reference model atmosphere with SHe = —1.0, es; = -3.75, ecr = -5.9, and ep e = -4.4. 
Lower plot: The emergent flux from the model atmospheres with modified abundance of individual elements minus the flux from a 
reference model. All fluxes were smoothed by a Gaussian filter with a dispersion of 10 A to show the changes in continuum, which 
are important for SED variability. The passbands of the uvby photometric system are also shown in the graph (grey areas). 
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3000 4000 5000 6000 
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Fig. 3. Magnitude difference Amj between the emergent fluxes 
calculated with an enhanced abundance of individual elements 
and the reference flux Hf 1 (see Eq.|6]l. The fluxes were smoothed 
with a Gaussian filter with a dispersion of 100 A. 



ence decreases with increasing wavelength. However, the be- 
haviour of the flux calculated for modified helium is different. 
The maxima of the relative brightness at the positions of the hy- 
drogen lines (especially close to the Balmer jump) remained un- 
detected in the previous analysis. These maxima are caused by 
the strengthening of the Lorentz wings of the hydrogen lines ow- 
ing to the higher density in the line-forming region in the models 
with higher helium abundance. It has not escaped our attention 
that the minimum of the relative magnitude difference at about 
5200 A caused by an accumulation of iron lines can be connected 
with the well-known flux depression at these wavelengths (as 
discu ssed already by, e.g.. lKhan & Shu lvak 2007; Krticka et al. 
2009). 





He 

Si 

Cr 

\ Fe 














-0.2 0.2 0.4 0.6 


0.8 1 1.2 



b 



phase tp 

Fig. 4. Predicted light variations of CU Vir in the Stromgren pho- 
tometric system calculated using abundance maps of one ele- 
ment only. The abundance of other elements was fixed. Light 
curves in individual filters were vertically shifted to better 
demonstrate the light variability. 



4. Predicted light variations 

Predicted light c urves are calcula t ed fro m the surface abundance 
maps derived by iKuschnig et al.l (1999) and from the emergent 
fluxes computed with the SYNSPEC code, applying Eq. [5] for 
individual rotational phases. 
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-0.02 



0.4 0.6 
phase <p 

Fig. 5. Predicted light variations of CU Vir (solid lines) com- 
puted taking into account the helium, silicon, chromium, and 
iron s urface abundance distributions derived by Kuschnig et al. 
(1 19991). The obser ved light variations (dots) are taken from 
Pyper et all dl998l) . The light curves in individual filters were 
shifted vertically to better demonstrate the light variability. 



To study the influence of individual elements separately, we 
first calculated the light variations with the abundance map of 
one element only (Fig. assuming a fixed abundance of other 
elements (e He = -1-0, esi = -3.75, S& = -5.9, sp e = -4.4). 
From Fig.|4]it follows that iron, silicon, and chromium contribute 
most to the light variations, while the contribution of helium is 
only marginal. This is because of the large overabundance of 
these elements in the spots and by their large abundance vari- 
ations on the stellar surface. The amplitude of the light varia- 
tions increases with decreasing wavelength, as can be expected 
from the plot of the magnitude difference Am^ in Fig. [3] Because 
the overabundant regions are brighter in the uvby colours, the 
predicted light variat ions reflect the equivalent width variations 
(iKuschnig et al ll999l Fig. 1). The light maximum occurs at the 
same phase at which the equivalent widths of a given element are 
the largest. Because this happens at slightly different phases for 
individual elements, the light curves in Fig.|4]are slightly shifted. 

Taking into account the surface distribution of helium, sil- 
icon, chromium, and iron in the calculation of the light curves 
(Fig. |5}, we obtained a good agreement between the observed 
and predicted light curves in the v, b and y bands of the 
Stromgren photometric system. On the other hand, our models 
are able to explain only about half of the amplitude in the u fil- 
ter. The disagreement between the predicted and observed light 
curves is mostly apparent around phase <p — 0.6. Note also that 
a similar disagreement visible in u can be also found in other 
filters, but to a much smaller extent. These differences clearly 
point to an existence of some additional, unknown mechanism 
working especially in the violet band that still needs to be inves- 
tigated (Fig. |5] and see also Sect. |7). The discrepancies between 
the predicted and observed light curves increase when compar- 
ing the predicted and observed colour indices (see Fig.|6]l. While 



0.02 



£ 0.04 



0.08 



0.12 





v-b 



-0.2 0.2 0.4 0.6 0.8 1 1.2 
phase <fi 

Fig. 6. Predicted variations of colour indices (solid lines) cal- 
culated from the helium, silicon, chromium, and iron surface 
abundance maps compared wit h the observa t ions. Observed light 
variations (dots) are taken from lPvper et al.l d 1998b . Light curves 
in individual filters were vertically shifted to better demonstrate 
the light variability. 

the (b—y) data agree reasonably well, the (v-b) curves are mutu- 
ally shifted, and the predicted metallic index nt\ = (v-b) — (b—y) 
shows a significantly lower amplitude than the observed one. 

The inhomogeneous surface distribution of individual ele- 
ments causes bright spots on the stellar surface. The spots, whose 
surface distribution can be derived using abundance maps and 
model atmospheres (see Fig. [7j>, cause the light variability. 

5. Ultraviolet variations 

We have shown that the light variability of CU Vir is caused by 
the redistribution of flux from the far UV to the near UV and 
visible regions. Consequently, the light variability in the far UV 
region should be in antiphase with the visual one. This behaviour 
was indee d found in a de tailed analysis of IUE observations of 
CU Vir bv lSokoiovl d2000b . 

To test these predictions quantitatively as well, we ex- 
tracted IUE observatio ns of CU Vir from the INES database 
dWamsteker et al.ll2000i see Table [ATTb using the SPLAT pack- 
age dDraperl 120041, see also ISkodall2008l) . Here we used low- 
dispersion large aperture spectra in the domains 1250-1900 A 
(SWP camera) and 2000-3000 A (LWR camera). 

5.1. Narrow-band UV variations 

As a first comparison of the UV fluxes we concentrated on 
narrow-band variations. For this purpose we smoothed the ob- 
served and predicted fluxes with a Gaussian filter with a disper- 
sion of 10 A. 
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Fig. 7. Emergent intensity from individual surface elements of CU Vir at various rotational phases. Upper panel: visible y band. 
Lower panel: UV band centred at 1250 A. Both for ji — 1 . 
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Fig. 8. Comparison of the predicted flux (mean and its variation, blue) with corresponding observed quantities (black). Both pre- 
dicted and observed fluxes (IUE) were smoothed by a Gaussian filter with a dispersion of 10 A. 



The resulting predicted and observed UV fluxes are given in 
Fig. [8] To avoid possible problems with absolute IUE calibra- 
tion, we normalised the predicted fluxes by a multiplicative fac- 
tor, which yields the best match between observations and pre- 
diction in Fig. [8] The factor was kept fixed for all wavelengths 
in all subsequent calculations. Generally, the mean predicted 
and observed fluxes agree well, with some minor differences. 
The mean observed flux is slightly lower than the predicted one 
in the region 1600 - 1850 A, whereas it is slightly higher in 
2100 - 3000 A. The amplitude of the observed flux variations 
agrees well with the observed one in the regions 1250 - 1450 A 
and 2550 - 3000 A. In the remaining regions the predicted am- 
plitude is lower than the observed one. 

As can be seen from Fig. [9] the predicted and observed flux 
variations in the selected wavelengths agree very well. Our sim- 
ulations are able to explain most of the observed features in the 
UV light curve. The amplitude is highest in the far-UV region 
1250 - 1400 A, reaching nearly 1 mag at 1250 A. The light vari- 
ation s in this region are mainly caused by silicon (c.f., Fig. [2] 
and Sokolov 120061 120101) . The overall agreement between the 
predicted and observed light curves in this region indicates that 
silicon abundances are well mapped. 

The amplitude of the light variations is very low in the re- 
gion between 1600 - 1900 A. The silicon-rich patches are dark 
in this wavelength region, whereas the iron-rich ones are bright, 



causing a near cancelation of any light variability in common. 
Note, however, that a fine structure of observed variations is not 
completely reproduced by the models, indicating either that the 
model atmospheres need to be improved or, which is even more 
likely, that the light variability has other sources. 

Our model s are able to nicely reproduce the observed 
dSokolov| [2000) antiphase variations in the far-UV on one side 
and the near-UV and visible region on the other side (see also 
Fig- ID- However, as was already clear from Fig. [8] the observed 
light amplitude is higher than the predicted one in the region 
2000 - 2500 A. This disagreement together with the difference 
in the u light curves (see Fig. [5]) indicates a presence of an ad- 
ditional now unidentified source of light variability. This might 
be connect ed with a chemical elem ent whose abundance was not 
mapped by iKuschnig et alJ (Q~999). This element, together with 
silicon, might contribute to the light variability in this region. 

Interestingly, despite the disagreement of the observed and 
predicted light variations in the u colour of the Stromgren pho- 
tometric system (Fig.|5]l, our models are able to nicely reproduce 
most of the light variations in the region 2550 - 3000 A (Pig. [9j>. 
The amplitudes of the observed light curves in this region and 
their shapes are a result of the inhomogeneous surface distri - 
bution of silicon, iron, and chromium (see also Sokolov 2010). 
Even more subtle effects, like the mutual shift of the light max- 
ima at wavelengths 2550 A and 2800 A, can be explained by our 
models. 
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Fig. 9. Comparison of the predicted (solid line) and observed (dots) UV light variations for different wavelengths. Curves for 
individual wavelengths were vertically shifted to better demonstrate the variability. 
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As can be seen from Fig. [10] the mean flux from CU Vir on wavelengths lower than 2400 A, whereas it is brighter than 
does not correspond to the flux calculated for the solar chemi- normal stars in the near-UV and visible regions, 
cal composition. The redistribution of the flux from the far-UV 
to near-UV and visible regions is apparent even when compar- 
ing the average fluxes. Consequently, CU Vir is fainter than the 5Z Monochromatic variations 

normal stars with the same effective temperature in the far-UV A i thoU g h ^ compar i SO n of the narrow-band variations revealed 

the regions where the disagreement between the observed and 
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predicted flux variations occurs, the narrow -band variations are sion of 1 .3 A, which roughly corresponds to a broadening of IUE 
inadequate to figure out the origin of these variations. For this data, and compared it with observed flux variations. The result- 
purpose the monochromatic variations are much more conve- ing monochromatic UV fluxes are given in Fig.fTTl 

nient. To compare the monochromatic fluxes, we smoothed the . f . ... ■ , r • c - rrn 

, _ F . . ., _ . ' ., ,. Most features that appear as individual lines in Fig. II Hare in 
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F F fact blends of a large number of individual lines (iron is a typi- 
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cal case). Only in a few cases it is possible to identify individual 
lines. The most numerous individual lines are those of silicon. 
The monochromatic fluxes in Fig.[TT]show that even though the 
continuum far-UV flux in the region A < 1500 A is relatively 
well fitted, there are some differences in the strengths of sili- 
con lines. That the strength of the Cn 1335 A doublet varies 
in antiphase with silicon indicates that silicon-rich regions are 
carbon-poor. From the weakness of the Alii 1671 A line we can 
conclude that this element is significantly underabundant with 
respect to the solar value. 

The predicted and observed fluxes during the light minimum 
close to the phase = agree very well. However, for phases 
(p » 0.5 the observed and predicted fluxes disagree in some re- 
gions, indicating that an additional opacity source is operating 
here, leading to yet another redistribution of the flux from far- 
UV to near-UV and visible regions. From Fig.QT|it is also pos- 
sible to identify these missing features that are responsible for 
the unexplained part of the light variability. From the plot for 
= 0.54 we can conclude that the missing opacity sources are 
located in the wavelength interval 1350 - 1800 A. 



6. Detailed analysis of observed and simulated light 
curves 

6. 1. Description of observed light curves 

The detailed analysis of the observed light curves in the wave- 
length region from 1250 A to 7600 A was made on basis of all 
available photometric data of CU Vir, inc luding very prec i se and 
relia ble data ob t ained in 1987-1997 by lAdelman et ail fl992) 
and [Pvper et al. (1998), and the data derived from IUE spec- 
trophotometry (see Sect.|5]l. The complete list of all photometric 
observations taken in 53 photometric bands that are quite evenly 
distributed along th e whole spectral i nterva l obtained in 1955- 
201 1 is published in Mikulas ek et alJ(l201 ll) . For the calculation 
of photometric phases we used a new ephemeris that takes into 
account the long-term variability of the period (see Sect. 12. U . 
The quantity and the quality (the mean weighted uncertainty 
of an individual photometric measurement is 5.3 mmag) of this 
photometric material enable us to investigate the photometric be- 
haviour of CU Vir with an unprecedented accuracy. 

First we studied the properties of the light curves in indi- 
vidual passbands. Our results show that the shapes of the light 
curves remained constant during the last half century. Each of 
them can be well expressed by a smooth single wave curve - a 
harmonic function of a low order. Nevertheless, the shapes of the 
light curves in different bands are apparently different. Applying 
the weighted advanced principle component analysis (APCA) to 
the parameters of the harmonic fits (which allows one to find hid- 
den relati onships among th em, for details see e.g. a brief intro- 
duction in Mikulasek 2007) we arrived at two principal conclu- 
sions: 1) Each light curve studied can be satisfactorily well fitted 
by a second-order harmonic polynomial - the amplitudes of the 
third and higher harmonics are always bellow 0.5 mmag. 2) All 
light curves studied can be well expressed by a linear combina- 
tion of only two basic light curves Fi((p), F2OP), the amplitude 
of the third and higher principal components do not exceed 1.0 
mmag. 

This allows us to build a relatively simple two-component 
phenomenological model with a minimum of free parameters 
valid for all studied light curves, where each of them is expressed 
by a linear combination of two principle functions Fi(<p), F 2(0) 




Fig. 12. Comparison of the observed effective semiamplitudes 
A(A) (see Eq.|71 dots with error bars) and semiamplitudes of sim- 
ulated light variations (small dots connected by solid line). 
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Fig. 13. Comparison of the form of observed light curves de- 
scribed by the dimensionless value if/ (see Eq. [7] dots with error 
bars) and with the value predicted by our models (small dots 
connected by solid line). 



of a rotational phase (f> 



Am(<p, A) = A(A) 



sin (^0W)) F i(0) + cos F 2 (4>) 



(7) 



where A(A) is the effective semiamplitude in the band cen- 
tred on the wavelength A, i[/{A) is the parameter explicitly 
determining the shape of the light curve in the wavelength 
A. The functions F\(<p), F2(<f>) are determined by parameters 
7t, 72, 0n, 012, 021, and 022 

F;(0) = cos(y/) cos[2?r(0-0n)] + sin(y,) cos[4 7r (0-0 Q )L (8) 

where I =1 and 2. Functions F\(<p), ^(0) were determined by 
means of APCA, where we confined ourselves to the first two 
principle components. They represent the pair of normalised mu- 
tually orthogonal vectors in the 4-D space of the Fourier coeffi- 
cients. Eq. [7] then means that the light curves defined by their 
vector in the space of the Fourier coefficients should lay in the 
plane determined by these two vectors of F\ (0), ^(0) functions. 

The light curves reach the maximum amplitude for wave- 
lengths shorter than 1600 A, whereas their shapes are nearly 
identical in this region (see Figs.|T2j[T3] an d Table lA.21 . That is 
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why we rotated the orthogonal vector base in its plane by an an- 
gle of -0. 12 rad so that the x-axis intersects these points. The pa- 
rameters of the principle light curves F\(<p), F2(<p) are then j\ = 
0.1244 rad, <f> u = -0.0132, 12 = 0.00442, y 2 = 0.2468 rad, 
021 = -0.2655, and <fe = 0.0912. 

6.2. Comparison of observed and predicted light curves 

The shapes of the simulated light curves can be described by 
the same procedure as the observed ones, consequently we can 
compare them unambiguously. We have arrived at the following 
conclusions: 

1. The shapes of both the observed and the simulated light 
curves strongly depend on the wavelength. There are differ- 
ences not only in the amplitudes of variations A(A) but also 
in their shapes, quantified by the parameter ifr(A). This can be 
understood as the result of the influence of several elements 
(at least Si, Fe, Cr). 

2. The shapes of the basic phase curves F\(<f>), F 2 (0) and the 
dependence of their parameters A(A) and tfr(A) on wavelength 
can be caused by two vast overlapping photometric spots 
centred at phases <p m = 0.3209(24) and <f> 02 = 0.6244(23) 
with different contrasts on the stellar surface. These spots 
can be identified as spots with overabundant Si and iron-peak 
elements (Cr, Fe). 

3. The uneven colours of the photometric spots described by 
the parameter tfr are illustrated in Fig. [9] with the IUE light 
curves in wavelengths 1500 A, 2100 A, and 3000 A. 

4. The dissimilar spectral energy distribution in spots is the rea- 
son why we do not find any 'null' point with zero amplitude 
of the flux variations. 

5. The simulation of the light variations cannot fully explain 
the observed amplitudes in the whole studied spectral region. 
The simulated amplitudes are on average smaller by a few 
tens of percent than the observed ones (see Fig.fT2l. which 
implies that we have neglected some of the important sources 
of the light variations in our modelling. 



a good agreement between predicted and observed SED vari- 
ations in the far-UV region, where the silicon dominates. On 
the other hand, the predicted and observed equivalen t widths of 
chromium and iron in Fig. 1 of Kuschni g et al](fl9 99) disagree to 
some extent. Our tests showed that especially the modified abun- 
dance of iron could help to explain some part of the discrepancy 
between theory and observations. 

Another source of the discrepancies might be connected with 
limitations of the model atmospheres used for the abundance 
analysis. All these considerations point to a need of new more 
precise abundance maps of CU Vir. 

7.2. Influence of additional elements 

We have shown that most of the light variability of CU Vir 
is caused by the uneven distribution of chemical elements. 
Consequently, it is possible that a part of the remaining dis- 
crepancy between theory and observation is also caused by 
ot her element(s), whose surface distribution was not mapped 
by Kuschnig et al. ( 1999). Both the optical and UV light curves 
provide strong constraints on the opacity caused by this as yet 
unidentified element, which might redistribute the flux from 
UV region of about 2000 — 2500 A, especially to the region of 
Stromgren u. 

We tested if either of the chemical elements currently in- 
cluded in the TLUSTY model atmospheres could cause these 
light variations. Excluding magnesium because of its low abun- 
dance previously, no other element included in TLUSTY (i.e., 
C, N, O, Ne, Al, and S) is able to cause the remaining light vari- 
ations observed in CU Vir. Consequently, it is likely that another 
element (especially the iron-peak ones) could be the cause. 

Very recent observations in a broad spectral range suggest 
that titanium and oxygen could also be contributors to inho- 
mogeneous abundance structures on the surface of CU Vir. 
Especially titanium (provided it is significantly overabundant) 
is one of the potential causes of the UV variations we cannot 
simulate so far. 



7. Discussion 



7.3. Vertical abundance stratification 



Despite the overall good agreement between the predicted and 
observed light curves, there still remain some differences be- 
tween the simulations and observations. These are especially ap- 
parent in the Stromgren u filter and in the UV region 2000 - 
2500 A. Note, however, that CU Vir has the highest amplitude of 
all stars modelled in detail so far, consequently it is encourag- 
ing that even this high amplitude can be explained for the most 
part. Here we discuss some possible reasons for the remaining 
disagreements. 

7.1. Limitations of abundance maps 

The limitations of abundance maps used could be an im- 
portant source of discrepancy between prediction and ob- 
servation. While the fine structure of abundance maps does 
not significantly i nfluen ce the predicted variability (as shown 
bv iRrticka et al] 12 009 using different abundance maps of 
iLehmann et alJ 20071) . other effects may be significant. The pre- 
dicted SED variations are sensitive to the maximum abundance 
and abundance amplitude of a given element in the map. The 
observed variations of silicon equivalent widths are nicely re- 
produ ced by the abundance maps (see Fig. 1 of iKuschnig et alJ 
1999), indicating their high reliability. This is also supported by 



Verti cal abundance strat i fication is obser ved in some CP stars 
(e.g., Ryabchikova 2004). ISokolovl d2010h proposed that the ver- 
tical abundance stratification may influence the light variability. 
Our test calculations confirmed these expectations. The models 
with overabundant iron in the outer regions (for the Rosseland 
optical depth tr oss < 0.1) indeed show a larger magnitude dif- 
ference in u than in the v, b, and y colours. Consequently, ver- 
tical abundance stratification could possibly explain the differ- 
ence between the phases of maxima in individual Stromgren fil- 
ters. However, the influence of vertical abundance stratification 
on the UV region is relatively low, consequently another pro- 
cess is needed to explain the difference between observation and 
theory in this region. 

7.4. Surface temperature variations 

Surface temperature differences and variable temperature gra- 
dients were also suggested as poss i ble causes of the CP star 
light variability (IWeiss et all Il976t IStepienl 119781). Hot star s 
may retain subsurface convection zones (ICantiello et alJ l2010). 
which can generate local magnetic fields. These magnetic fields 
may give rise to the surface temperature differences, and con- 
sequently cause the light variability. However, mild differences 
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between the observed and predicted light curves do not indicate 
any (effective) temperature differences on the surface of CU Vir. 
This is supported also by the fact that e. g. the predicted mean 
of the flux distribution simulates the observed one fairly well, 
as depicted in Fig. [8] Moreover, a good agreement between the 
observed and predicted light curves in the far-UV region and the 
vby light curves in the visible region excludes the temperature 
differences as a main source of the light variability. 

7.5. influence of the turbulent velocity 

In our study we assumed a generic value of the microturbu- 
lent velocity 2kms~'. This parameter, which roughly accounts 
for atmospheric velocity fi elds, likely has a zero value in cor- 
responding normal stars (Lands treet et al.l 120091) . We kept a 
nonzero value here as a very rough approximation of Zeeman 
line broadening. The adopted value of the microturbulent veloc- 
ity may, however, influence the emergent flux. At higher micro- 
turbulent velocities the line transitions are able to absorb radia- 
tion more effectively, increasing thus the temperature in the con- 
tinuum forming regions. To estimate the magnitude of this effect, 
we calculated an additional model with higher microturbulent 
velocity 4kms~ I and assuming enhanced abundance of heavier 
elements (su e = —1.0, esi = -2.25, sq- = -4.9, ep e = -3.4). We 
compared the resulting flux distribution with the model with the 
same chemical composition, but with a standard microturbulent 
velocity of 2 km s _1 . The calculated magnitude difference (Eq.O 
between these models has its minimum -0.05 mag in the near- 
UV region 3000 - 3800 A and a maximum about -0.10 mag in 
the region 2200 - 2550 A. These are the regions where the most 
apparent differences between observed and predicted light vari- 
ations occur. Consequently, a higher value of the microturbulent 
velocity and/or surface microturbulent velocity distribution (see 
Sect. 17. 8t cannot be ruled out as a possible cause of the remain- 
ing difference between theory and observation. 

7.6. The effects of fast rotation 

CU Vir belongs to the fast rotators among CP stars. From TableQ] 
we can infer its rotational velocity v rot = 320 km s _I , indicating a 
rotational velocity close to the critical one. To quantify this, the 
stellar radius has to be known with sufficiently high precision. 
Using the stellar parameters derived from spectroscopy and pho- 
tometry (i.e., v lot sin/, inclination, and period, see Table [TJ, we 
can estimate the equatorial radius to be ,R eq = 3.3 + 0.6 R e . This 
value agrees well with the stellar radius d erived from the evolu - 
tionary tracks in the T e ff - \ogg plane of ISchaller et alJ ([1992), 
which isR = 3.2+0.5 R Q (the derived mass isM = 3.8 ±0.2 M ). 
Note, however, that this yields a significantly higher radius than 
that derived from the observed UV flux, which is R = 1.9+0.1/?© 
assuming a distance of 79+1 pc (I van Leeuwenl2007l) . This possi- 
bly points either to a problem with the absolute flux calibration, 
or to an inclination that is too low. The latter is supported by a 
lower radius of R = 2 .3 + 0.1 R n derived from photo metry and 
evolutionary tracks by Kochukhov & Bagnulol (l2006h . 

To study the effect of fast rotation at its extremum, we as- 
sumed the equatorial radius R eq = 3.3 R e , and M = 3.8 M , 

the polar radius is then R p = R eq (l + v^ ot R eq / (2GM)) = 
2.1 R Q dCollinsI [19631) . The ratio of the rotational velocity to 
the critical one is then v rot /vimt = v mil ^2GM/3R p = 0.75. 
If the star rotates this rapidly, the polar to equator difference 
in the effective surface gravity and effective temperature are 



A logs = log {GM/Rl) - log(GM/^-v^ t /geqj « -0 -4 and 
Ar eff » 3000 K (assuming 7W(ff) ~ g 1/4 (fl). lvon . Zeipellll924t 
ISlettebakl[T949l: lOwocki etal.ll 1 9941) . 

The variations of local surface gravity and effective tempera- 
ture connected with a fast rotation by itself cannot raise any light 
variability (assuming fixed axis of rotation) because of their axial 
symmetry. However, even axisymmetric surface variations may 
modify the light variations caused by inhomogeneous elemental 
surface distribution via several effects. First, surface layers with 
different abundances may respond differently to the temperature 
and surface gravity variations, modifying the flux distribution. 
Second, the area of the surface element is modified through the 
oblateness of the surface of a rotating star. Finally, the direction 
of the beam pointing to the observer relative to the local outward 
normal and its cross-section are different on spherical and oblate 
surface. 

To test the influence of these effects, we calculated an ad- 
ditional grid of model atmospheres corresponding to the stel- 
lar equator {T e ff fiq = 11 500 K, logg eq = 3.71) and to the pole 
(Teff.p = 14 300 K, loggp = 4.16) with different abundances of 
silicon and iron (after Table 0. The polar and equatorial effec- 
tive temperature and gravity were c alculated assum i ng th at the 
CU Vir parameters r e ff, g derived by Kuschnig et al. ( 1999) rep- 
resent some kind of mean over the stellar surface and that these 
parameters correspond to some particular surface region. In this 
case T e g(&) = [gity/g] 1 ^ T e g. For simplicity we assumed a fixed 
abundance of helium and chromium here (en e = — 1, £<> = 5.9) 
because these elements are not the main sources of the light vari- 
ability. The resulting light curve was derived using Eq. [5] now 
interpolating the intensities in Eq. [4] also between the models 
corresponding to the pole and to the equator. 

The inclusion of the gravity darke ning only (with the 
latitu de-dependent radius calculated after Harrington & Collinsl 
1968) leads to the decrease of predicted amplitude of light varia- 
tions by up to 0.01 mag. This is because for both silicon and iron 
the values of Am^ Eq. [5]decrease with decreasing temperature. 

The inclusion of different surface areas caused by the oblate- 
ness of the star leads to another small modification of the light 
curve, which is mostly negligible (the difference is up to one 
millimagnitude). We did not consider the effect of the difference 
of beams pointing to the observer (in an oblate star compared to 
the spherical one), but given the small influence of other effects, 
we expect that this is also negligible. 

We conclude that the effect of the gravity darkening modi- 
fies the light curves, but it does not seem to be the main reason 
for the discrepancy between predicted and observed light curves. 
But the effect of gravity darkening on the light curves introduces 
a new possibility to test the theory of gravity darkening in the 
future. To this end, the influence of this effect on the abundances 
derived from the Doppler imaging has also to be accounted for. 

7. 7. Influence of the magnetic field 

CU Vir is known to host a large-scale magnetic field, which 
could, in principle, contribute to the light variability via the rota- 
tional modulation of the magnetic field intensity and thus opacity 
in Zeema n broadened spectral features. Assuming a dipolar field 
geometry Trig ilio et al.l (2000) based on p hase-resolved longi- 
tudina l magnetic field measurements from IBorra & Lands treet 
(11980 ) determined a polar magnetic field B p = 3 kG and the 
angle between the stellar rotation axis and the line-of-sight 
/3 = 74°. Regarding magnetic CP stars, this is a fairly moderate 
magnetic field and, as demonstrated bv lKhan & Shulvakl {2006) 
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based on detailed model atmosphere calculations with anoma- 
lous Zeeman effect and polarized radiative transfer, such a field 
does not produce significant effects on parameters observed in 
different photometric systems. In addition, with the inclination 
angle i = 30° the surface average magnetic field modulus varies 
only by A ( \B\ ) « 300 G during the rotational cycle, which is far 
too small to induce any detectable changes in light-curves. 

On the other hand, as shown by Shul vak et alj d2008l) . the 
combined impact of the magnetic field and the realistic chem- 
istry could be more important than the effect of using only in- 
dividual abundances. Thus, the effect of including the magnetic 
field in the computations of photometric parameters would be 
more stronger in regions of spots with enhanced abundances 
compared to the case of a homogeneously stratified atmosphere. 

Because we are mostly interested in the estimate of the max- 
imum possible amplitude that the magnetic field could introduce 
in the light curves relative to the inhomogeneously distributed 
abundances, we computed se veral magnetic mod el atmospheres 
using the LLmodels code (Shulvak et al. 2004). Anomalous 
Zeeman splitti ng and polarised radiativ e transfer were included 
as described in lKhan & Shulvakl (12006). We found that the high- 
est impact of the magnetic field for a spot with enhanced Fe, 
Si, and Cr, located at the pole with \B\ = 3 kG and equa- 
tor with \B\ = 1.5 kG amounts to Am « 0.005 mag. This 
difference subsequently decreases for other Stromgren filters. 
Accounting for the surface averaged magnetic field modulus re- 
sults in A ( \B\ ) = 2.2 kG at the phase when the magnetic pole 
is visible and A{ \B\ ) = 1.9 kG at the phase of magnetic equator 
respectively. The effect is then reduced to Am » 0.001 mag. All 
this is of an order of magnitude less than needed to explain the 
deviations between observed and predicted non-magnetic ampli- 
tudes of the M-parameter shown in Fig. [5] We therefore conclude 
that the magnetic field has little or negligible effect on the light- 
curve appearance and can be ignored in the present study. 

7.8. Convection zone in helium-rich models 

In addition to the central convection zone, a hot star may have 
subsurface iron and helium convective zones dMaederl [l980: 
iMaederet al.l2 008: Cant iello et~alll2010h . Our models show that 
for helium-rich models with su e ^ 0, the helium convection 
zone moves towards the stellar surface and its top is located 
in the model atmosphere for Rosseland optical depths tr oss w 
1 - 100. The existence of a subsurface convecti on zone may 
have i nteresting astrophysic al consequences (e.g.. lMaeder et al] 
2008; Cantiell o et alJ l20loh . The convective zone may create 
a dynamo, generating the chromospheric activity and conse- 
quently also X- ray emission, whose existe nce is still puzzling 
in A type stars ( Schroder & Schmitt 2007). Moreover, convec- 
tion may cause surface turbulence, leading to inhomogeneous 
distribution of the turbulent velocity provided helium is also dis- 
tributed inhomogeneously. 

However, these effects are likely strongl y damped in CU Vir, 
which has a strong surface magnetic field (Landstre et & Borral 
119771) . Our results show that the magnetic field energy density 
dominates over the gas energy density to the optical depths of 
about tr oss * 10 3 , suppressing any motion perpendicular to the 
magnetic field lines. 

7.9. NLTE effects 

Although the TLUSTY model atmosphere code enables us to 
calculate NLTE models, we confined ourselves to LTE models, 



because we expect NLTE effects to be marginal for the light 
variability. To test this we calculated additional NLTE models. 
Because we do not have sufficient atomic data to calculate NLTE 
models including chromium, we forced LTE for this element, 
even if for the remaining elements we assumed NLTE. 

We calculated an NLTE model with enhanced abundance 
of heavier elements (sue = -1.0, esi = -2.25, £& = -4.9, 
£Fe = -3.4) and compared it with the corresponding LTE model. 
The resulting fluxes differ by about 1 - 2 % in the visible and 
near-UV regions. The most pronounced changes appear in the 
far-UV region with A < 1500 A. Consequently, the changes ow- 
ing to NLTE are significantly lower than those owing to variable 
abundances, and cannot be the main source of the difference be- 
tween theory and observations. 

The NLTE effects are, however, significant in lines. The 
NLTE effects lead not only to the well-known strengthe ning 
of the core of hydrogen lines (e.g.. lLanz & Hubenv| [2007). but 
the line equivalent widths of other elements are also slightly af- 
fected. Finally, in NLTE some infrared lines appear in emission. 

8. Conclusions 

We successfully simulated the UV and visual SED variability of 
the helium-weak star CU Vir. We assumed that the light vari- 
ability is caused by the inhomogeneous surface distribution of 
elements and used model atmospheres to predict the light vari- 
ability. 

Individual chemical elements are distribu ted inhomoge- 
neousl y on the surface of CU Vir, as derived bv lKuschnig et al.l 
(1 19991) by Doppler mapping. The chemical composition influ- 
ences the emergent flux through flux redistribution from the far- 
UV to near-UV and visible regions. The bound-free transitions 
of silicon and the bound-bound transitions of iron and chromium 
mostly cause the flux redistribution in CU Vir. As a result of the 
flux redistribution, the individual surface elements display dif- 
ferent brightness in individual photometric bands, although the 
total (frequency integrated) emergent flux is the same for all sur- 
face elements. The inhomogeneous surface brightness manifests 
itself by the light variations caused by stellar rotation. 

The inhomogeneous surface distribution of silicon, 
chromium, and iron is able to explain most of the observed 
UV and visible SED variations. We successfully reproduced 
the antiphase behaviour of the light curves in the far-UV and 
visible regions. We emphasise that the variability seen in 
the visible is just a faint gleam of the variability seen in the 
UV. While the amplitude of the light curves merely reaches 
about a few hundredths of magnitude in the visual domain, it 
reaches about 1 mag in the UV. However, our models are able 
to reproduce just part of the variability seen in the UV region 
2000 - 2500 A by IUE and in the u filter of the Stromgren 
photometric system. Another mechanism(s) has to be invoked to 
explain this difference between the observed and predicted light 
curves. This so far unidentified mechanism could be connected 
with inhomogeneous distribution of an element whose surface 
abund ance distribution was not mapped by Kuschnig et al.l 
(1 19991) . but also other effects may contribute to this difference. 

Our models nicely reproduce the observed SED for the phase 
<f> ~ 0, when the regions with lowest elemental abundances are 
seen. On the other hand, there are some discrepancies between 
simulated and observed SED for the phase (p w 0.5, when the 
overabundant regions appear. These discrepancies have likely 
the same origin as the differences of the light curves. 

Our study provides additional evidence that the light vari- 
ability of chemically peculiar stars is mostly caused by the inho- 
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mogeneous surface distribution of individual chemical elements, 
flux redistribution, and stellar rotation. The comparison of ob- 
served and predicted SED and its variation may serve as a test of 
opacity sources included in current model atmospheres. Finally, 
as a byproduct, it provides an image of the stellar surface, which 
(besides the interferometry) is the only way how these images 
can be derived. 
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Appendix A: Long tables 

Table A.l. List of the IUE observations of CU Vir 



V .- Cll 1 1V_.± CI 
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2 400 000+ 


Phase 
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znSRzl 07 'SOS 


4^6 


LWR 


3448 


tJOOt^. 1 U OO 
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4^040 SR9R4 
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T WR 

l^VV IS. 


4047 






LWR 
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1.000 


SWP 


3864 


43883.91128 
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SWP 
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43883.95163 
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43883.99343 


0.300 


SWP 


3867 


43884.03658 


0.382 


SWP 


3868 


43884.07952 


0.465 
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43884.12081 


0.544 


SWP 


4670 


43949.80501 


0.695 


SWP 


4671 


43949.84709 


0.776 


SWP 


4672 


43949.88714 


0.853 


SWP 


4673 


43949.92677 


0.929 


SWP 


4674 


43949.96918 


0.010 



Table A.2. Observed and simulated effective semiamplitudes of 
light curves and parameters describing their forms in individual 
photometric bands of ultraviolet and optical regions (see Eq.[7j. 



band 
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Notes. Af denotes number of observations 
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